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Recent studies show that in interdependent networks a very small failure in 
one network may lead to catastrophic consequences. Above a critical fraction of 
interdependent nodes, even a single node failure can invoke cascading failures 
that may abruptly fragment the system, while below this "critical dependency" 
(CD) a failure of few nodes leads only to small damage to the system. So 
far, the research has been focused on interdependent random networks without 
space limitations. However, many real systems, such as power grids and the 
Internet, are not random but are spatially embedded. Here we analytically 
and numerically analyze the stability of systems consisting of interdependent 
spatially embedded networks modeled as lattice networks. Surprisingly, we find 
that in lattice systems, in contrast to non-embedded systems, there is no CD 
and any small fraction of interdependent nodes leads to an abrupt collapse. 
We show that this extreme vulnerability of very weakly coupled lattices is a 
consequence of the critical exponent describing the percolation transition of a 
single lattice. Our results are important for understanding the vulnerabilities 
and for designing robust interdependent spatial embedded networks. 

Complex systems, usually represented as complex networks, are rarely isolated but usually 
interdependent and interact with other systems l|-|3|]. Recently it was shown that a coupled 
networks system is considerably more vulnerable than its isolated component networks 



22] . A failure of nodes in one network leads to a failure of dependent nodes in other networks, 



which in turn may cause further damage to the first network, leading to cascading failures 
and catastrophic consequences. It was shown that the strength of the coupling between the 
networks, represented by the fraction q of interdependent nodes, determines the way the 
system collapses {5, Is, U^. For strong coupling, that is for high fraction of interdependent 
nodes, an initial damage can lead to cascading failures that yield an abrupt collapse of the 
system, in a form of a first-order phase transition. Reducing the coupling strength below a 
critical value, q c , leads to a change from an abrupt collapse to continuous decrease of the 
size of the network, in a form of a second-order phase transition. This new paradigm is in 
marked contrast to the common knowledge represented by a single network behavior. In a 
single network a failure of few nodes can make only a small damage to the network, while 
a system of interdependent networks might be functioning well with most nodes connected, 
but a further failure of even a single node may lead to a complete collapse of the entire 
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system. Thus, the existence of an abrupt collapse phenomena in interdependent networks 
makes such systems extremely risky. Thus, understanding this phenomena is critical for 
evaluating the systems' vulnerability and for designing robust infrastructures. 

Current models focus on interdependent networks where space restrictions are not con- 
sidered. Indeed, in some complex systems the spatial location of the nodes is not relevant or 



not even defined, such as in proteins interaction networks 2J-[26| and the World Wide Web 



3. 



However, in many real-world systems, such as power grid networks and computer 
networks, nodes and links are located in Euclidian two-dimensional space. The dimension 
and universality class of embedded networks, whose links have a characteristic length, is the 
same for all embedded networks j^j. For example, in power grid networks the links have 



30] . Thus, to 



a characteristic length since their lengths follow an exponential distribution 
obtain the main features of the system under failures, we model such a spatial embedding 
network two dimensional lattice. 

Here, we study the case of interdependent lattice networks where a fraction q of nodes 
in each lattice randomly depends on nodes in other networks. We find that in the case of 
coupled lattices q c = and any coupling q > leads to a first-order transition. We show that 
the origin for this extreme vulnerability of coupled lattices compared to random networks 
(q c > 0) lies in the critical behavior of percolation in a single lattice, which is characterized 



by the critical exponent /3 < 1, in contrast to random networks where {3 = 1 [3jj. Our 
theoretical and numerical approaches predict that a real-world system of interdependent 
spatially embedded networks which are characterized by < 1 will, for any q > 0, abruptly 
disintegrate. Since for percolation of lattice networks it is known that for any dimension 
d < 6, f3 < 1 [29[], we expect that also interdependent systems embedded in d = 3 (or any 
d < 6) will collapse abruptly for any finite q. 

In addition of studying a pair of interdependent lattices we also analyze the stability of a 
more general case of a starlike network-of-lattices, where each of the n — 1 peripheral lattices 
is coupled to the root lattice with the same coupling strength q. We find that the root lattice 
abruptly collapses (first-order transition) for any q and any n > 1, however, the percolation 
behavior of the peripheral lattices is somewhat counterintuitive: for small n the peripheral 
lattices also abruptly collapse together with the root lattice, while for large n the peripheral 
lattices may remain functional even after the root collapses. 

Our analytical considerations consist of four stages. First, we analyze the conditions for 
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FIG. 1: Abrupt collapse of interdependent lattices, a, The size of the giant component 
Poo of two interdependent lattice networks, each of size 4000 x 4000, at steady state after random 
failure of 1 —p of their nodes. The solid lines are the solutions of Eq. ([1]) and the symbols represent 
simulation results. In the case of q = (circles), no coupling between the networks, the transition is 
the conventional second-order percolation, while for any q > (triangle downs, q = 0.2; diamonds, 
q = 0.4; squares, q = 0.6; stars, q = 0.8; triangles up, q = 1) the collapse is abrupt in the form 
of first-order transition, b, The number of iterations (NOI) vs. p, triggered by removing at each 
step a fraction Ap = 10~ 3 of the nodes. In the case of first-order transition (abrupt collapse due 
to cascading of failures), the number of cascading failures, NOI, diverges for p —>Pc, providing a 
useful approach for detecting accurately the critical threshold pc in simulations [191 ]. c, The solid 
line is the solution of Eqs. ([1]) and ([2]) for pc versus q and the symbols represent results from 
simulations. Note that for the case of no coupling between the networks, q = 0, we obtain the 
known result of a single lattice »e = 0.593 and for the case of complete coupling q = 1 we obtain 



Pc = 0.826, in agreement with 



32]. 



having an abrupt transition (first-order percolation transition) in a system of any coupled 
networks system. Second, we find the critical coupling strength q c below which the transition 
becomes continuous (second-order). Third, we show that for the case of two coupled lattices 
q c = 0, thus the transition is of first-order for any q > 0. Finally, we generalize our results 
to the case of network-of-lattices. 
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I. SYMMETRIC INTERDEPENDENT NETWORKS 



We analyze here, for simplicity and without loss of generality, the symmetric case where 
the two networks, A and B, have the same degree distribution. A fraction q of the nodes 
of each network are dependent on nodes randomly selected in the other network under the 
assumption of no-feedback condition, meaning that if node A{ depends on node Bj , and Bj 
depends on , then k — i Q|. The analysis of the non-symmetric case and the case with 
feedback which yield similar results are discussed in the sections [III and III II respectively. In 
the symmetric case a fraction 1 — p of nodes of both networks is initially randomly removed. 
As a result, a certain fraction of nodes become disconnected from each network, while a 
fraction g(p) of the nodes in each network remains in the giant component. Each node that 
has been removed or disconnected from the giant component causes its dependent node in 
the other network to also fail. This leads to further disconnections in the other network and 
to cascading failures. The size of the networks' giant components at steady state is given 
by -Poo = xg(x), where x is the solution of the self consistent equation [§] 

x - p(l - q) = p 2 qg(x). (1) 

The function g(x) represents the probability of a node to be connected to the giant com- 
ponent after random removal of a fraction 1 — x of the nodes and can be obtained either 
analytically or numerically from the percolation behavior of a single network. In general, 
g(x) has a critical value at x = p c such that g(x < p c ) = while g(x > p c ) > and 
monotonically increases with x. For any given value of p, a graphical solution of Eq. ([I]) is 
given as the intersection of the straight line y = x — p(l — q) and the curve y = p 2 qg(x). 
The value of p where the line and the curve tangentially intersect, p = p^, corresponds to 
a discontinuity in the solution of x (see Fig. [2^l) . This leads also to a discontinuity in the 
solution of Poo{p) which abruptly jumps to zero as p slightly decreases. At any p < p^, x is 
smaller than the critical value p c (of a single network) and therefore Poo(p < p 1 ^) = 0. 

In Fig. [T^a) we show Poo{p) from theory, Eq. ([1]), and simulations as a function of p for 
various values of q. For q = the lattices are not coupled and the percolation transition 
is the conventional continuous second-order, while for q > the transition is discontinuous 
first-order. As seen, the theory and simulations are in excellent agreement. In order to find 
in simulations p^ accurately we evaluate the number of cascade iterations, (NOI) 



13, for 



each p, which has a sharp peak (diverges for infinite systems) at p^ as seen in Fig. []Jb). 



The condition for the abrupt jump is that additional to ([I]), the derivatives of both sides 
of Eq. with respect to x are equal, 

l=p 2 qg'(x). (2) 

Solving Eq. (pQ) together with Eq. (|2j) for a given q provides the abrupt jump transition point 
p% and which yields the size of the giant component just before collapsing, P^ = x^g{x^). 

The size of the giant component at criticality P^ depends on the coupling strength q 
such that reducing q leads to smaller value of x^ and thus smaller P£. As long as x M > p c , 
P^ > and -Poo(p) discontinuously jumps to zero at p just below p%. However, for a certain 
q, x M — >■ p c and the size of the jump becomes zero since P*^ — > (see Fig. [2b, d). In this 
case the percolation transition becomes continuous. Therefore, the critical dependency q c 
below which the discontinuous transition becomes continuous, must satisfies both Eq. ([I]) 
and Eq.flSJ) for x — Y p c , and is given by 

p(l - q c ) = p c (3) 

p 2 qcg'(Pc) = l. 

The set of equations ([3]) is a general condition for finding the transition point, q c , from 
discontinuous to continuous percolation transition in a symmetric case of any two networks. 
Interestingly, the transition point, q c , depends only on the single network behavior near 
criticality, g(x — p c ), and not on the entire shape of g(x). As seen in the following, in the 
case of random networks g'(p c ) is finite yielding a finite solution for q c . However, for the 
case of lattice networks the derivative of g(x) diverges at the critical point, g'(p c ) = oo, 
yielding q c = 0. Therefore, from Eq. flS} follows that any coupling between lattices leads to 
an abrupt first order transition. 

The behavior of the percolation order parameter of a single network near the critical 
point is defined by the critical exponent j3, where g(x) x ^. Pc = A(x — Pr) 13 - The divergence of 



g'(x) of 2d single lattice when x — > p c is since f3 = 5/36 < 1 2j], 33|, [3^]. In contrast, for ER 
networks (3 = 1 which yields finite value of g'(p c )- Thus, q c can be explicitly calculated for 
coupled ER networks with average degree (k) from Eqs. ([3]) based solely on the behavior at 
criticality, p c = 1/ (k) and A = 2{k), yielding q c = 1 — 4r(y2(fc) + 1 — 1). 

Figure [2] demonstrates schematically the possible solutions of Eq. <^ at criticality, p = 
p 1 *, where the value of x M yields the size of the system just before the abrupt collapse, 




Pc P Pc P 

FIG. 2: Schematic solution of Eq. ([!]) for p = p^ in the case of two coupled (a) lattices and (b) 

I I 

ER networks. The mixed ER-Lattice case is studied by Landau et. al. [35]. The left-hand-side 
(straight line) and right-hand-side (curve) of Eq. (pQ) are plotted versus x, where the solution 
at criticality (x = x^) is marked as solid circle. For the case of (a) two coupled lattices, since 
d'(Pc) °o the value of x^ is always larger than p c , thus, the transition is of first-order for any 
non-zero q value (c). However, for the case of (b) such as for two ER networks, g'(x) is finite for any 
x. Thus, for q < q c there is no tangential touching point between the straight line and the curve 
and the solution x continuously approaching p c as p decreases so the transition is of second-order 
(d). 

-^oo = x ^g{ x ^)- While for (b) coupled ER networks and q = q c the solution is x M = p c 
and P£ = 0, meaning that the system continuously disintegrates, for (a) coupled lattices, 
x M > p c and > for any q, thus, the system of coupled lattices always undergoes a first- 
order transition. Abrupt collapses are extremely risky since most of the system's elements 
just before the collapse may be functioning well and no sign of warning appear. This is in 
mark contrast to single network behavior where a failure of small fraction of nodes can lead 
only to a small damage to the network. Next, we study the size of the giant component of 
coupled lattices at criticality Poo(j>c) f° r different coupling strengthes q, using Eqs. (0Q) and 
([2]) as well as numerical simulations. As shown in Fig. [3j in the case of a coupled lattices 
system for any q > 0, Poo(Pc) > 0> that is the size of the collapse, is larger than zero for 
any coupling between the lattices. We also compare, in Fig. [3], the size of the collapse in 
a system of coupled lattices and in a system of coupled random ER networks with average 
degree (k) = 4, same as the degree of the lattice. Figure [3] demonstrates the significantly 
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FIG. 3: The size of the collapse in coupled lattices compared to coupled random net- 
works. Comparison of the size of the giant component at criticality -Poo(Pc) for coupled lattices 
(circles) and coupled random (ER) networks with average degree (k) = 4 (squares) versus the 
coupling strength q. While for random networks with q < q c = 0.5 the size of the networks at 
criticality is zero, in coupled lattice network the networks abruptly collapse for any finite q > 0. 
Note also the huge differences in the network sizes at the collapse transition. The coupled lattices 
collapse at a significantly larger giant components compared to the ER case. The solid line repre- 
sent the theory for coupled lattices given by Eqs. (pQ) and ([2]), and the circles are from simulations. 
The dashed line represents the theory for ER coupled networks derived in Parshani et. al. 

increased vulnerability of the coupled lattices system compared with the random networks 
system. For example in the coupled random networks, for q = q c = 0.5 the size of the system 
at criticality is zero, while in coupled lattices the collapse occurs when the giant component 
is about 1 /2 of the original network. 

Following a similar approach, we analyze in section [Til the non-symmetric case of two 
interdependent networks, where a fraction q^ of the nodes of network j depend on nodes 
of network %. In section |TTJ we show analytically that the first lattice that breaks down 
undergoes always a first-order percolation transition, for any finite values of q\2 and qi\. 
Yet, the percolation behavior of the second lattice, does depend on g 12 and falls into one of 
two scenarios: the second lattice either (i) abruptly collapses together with the first lattice 
at p^2 = Pci or (ii) undergoes a second-order transition at p^ 2 < Pci, but still there is an 
abrupt fall at p^i- 



S 



When the first lattice completely collapses at p = p^i the cumulative failures of the second 
lattice can be evaluated by an effective occupation probability 



Peff = P\ 




(4) 



When qi2 > 1 — ^f- then p e ff < p c and the two lattices will abruptly collapse together at 

Pel 

p = Pci, while for q 12 < 1 — ph, p e ff > Pc and the second lattice undergoes a second order 
transition. In the second scenario from Eq. (jlj) we obtain that 



The general results obtained analytically for two coupled lattices are also useful for ana- 
lyzing network-of-lattices. For any given structure of the network-of-lattices, the first lattice 
that breaks down will undergo a first-order-transition. The reason is that the behavior of 
all other lattices connected to this lattice can be reduced to a single "effective" percolation 
function g e ff(x) and the considerations of two coupled lattices are valid. Once one or sev- 
eral lattices collapse the survived lattices that become isolated will undergo a second-order 
transition. 

Next, we analyze analytically and numerically a system of starlike network-of-lattices. 
The system consists of n lattices with equal dependency, = q for i = 1 or j = 1 and 
qi 7 j = elsewhere. A fraction of 1 —p of the nodes is randomly removed from all the lattices, 
Pi = p. The steady state of the system is given by 



where the giant component of the root lattice is Pno,i = xig(xi) and the giant component 
of the n — 1 peripheral lattices is Poo,2 — %29(%2) In the case of n = 2 Eqs. © and (J7J) 
reduce to Eqs. ([I]) and ([3]). For n > 2, the root lattice will collapse first since it depends 
on all the other lattices, consequently, it undergoes first-order transition for any value of q. 
Figure H] shows the percolation behavior of the peripheral lattices in the starlike network- 
of-lattices. We find that for small n the peripheral lattices collapses together with the root 
lattice at p = p^ v However, for large n, the collapse of the root is followed by a discontinuous 
fall in Poo t 2 at p — p^ which then undergoes a second-order transition at p% 2 < p^. 



Pc2 = Pc/{l-Ql2)- 



(5) 



X\ = p[qpg(x 2 ) - q + l] n 

n-2 

x 2 = qp^Xi' 1 g(xi) +p(l - q) 



(6) 
(7) 
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FIG. 4: Percolation in a starlike network-of-lattices. a, The size of the peripheral lattices 
Poo2(p) in a starlike network of n lattices with q = 0.2 as obtained analytically from Eqs. © 
and ([7]) (solid lines) and from numerical simulations (symbols). The discontinuous jumps in the 
peripheral lattices occur when the root lattice collapses at p = p^. For small n, the peripheral 
lattices collapse together with the root, p^ 2 = P^n while for large n, the peripheral lattices remain 
functional even after the root lattice collapses and continuously disintegrate at p^ 2 < p^, which 
is given by Eq. ([5]). b, The critical threshold of the peripheral lattices p^ 2 versus the number of 
lattices n for different values of coupling strength q from bottom to top: q = 0.1, 1.5, 2, 2.5, 3. 
The dashed lines are the theory and the points are from simulations. 

II. NON-SYMMETRIC INTERDEPENDENT NETWORKS 

Here we analyze the mutual percolation of any two interdependent networks for the non- 
symmetric case and find the conditions for discontinuous behavior. The steady state is given 
by the set of two equations |8|] 

xi = PiP2q2ig2(x 2 ) +pi(l - <?2l) (8) 

X 2 = P2Piqi29l(Xl) +P2(1 - 912) (9) 

where the size of the giant component of network i is given by Poo,i = Xigi{xi), Pi is the 
initial random occupation of network i and is the fraction of nodes in network j that 
depend on network i. For given gi(x) and g2{x), Eqs. (jHJ) and (jHJ) can be solved for x% and 
X2, for any values of pi and g^ . We assume, for simplicity, that in both networks a fraction 
1 — p is initially randomly removed, p\ = p 2 = p. Still, since g i2 7^ c/21 an d gi 7^ g 2 the two 
networks generally will not disintegrate together. 
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We define network 1 as the network that first collapses when p = p%, thus, the solution 
of x\ discontinuously jumps from x\ > x c \ for p > p% to X\ < x c \ for p < p%. For p = p% 
x\ = x^ L . g 2 (x) can be expanded using Tailor series around x 2 , the value of x 2 at p — )■ p 1 ^, 

92(x) = g 2 (x 2 ) + g 2 (x 2 )(x -x 2 ) + ... (10) 

For p ->■ p^, Eq. ^ becomes, x x = p 2 q 2 \[g 2 (x 2 ) + g' 2 (x 2 )(x 2 - x 2 )\ + p(l - <?2i), or 

x 2 = mx\ + n (11) 

where m = -2 — , and n = x 2 — p(l — q 2 \) — ^tt?4. 
Equation ([9]) becomes 

ma ; 1 +n = p 2 gi2fi'i(a;i)+p(l-gi2). (12) 

At the critical point, where p = p^ the solution x\ = x^ satisfies both Eq. (|T2|) and the 
tangential condition, obtained by taking the derivative of Eq. (fT2l) with respect to x\ 

m = p 2 q 12 g[(x 1 ), (13) 

which then becomes 

p^q^ig'iix^g'^) = 1. (14) 

The solution of Eqs. f|T2|) and ffT4j) provides p(f and a;^ of the critical point of the first-order 
transition. Decreasing the values of q\ 2 and 521 leads to lower value for x± and, consequently, 
lower value for Poo,i{Pc) = x i9i( x i)> the size of the giant component at criticality. If 
Xi — > x c i then Poo,i(Pc) — an d the transition becomes second-order. 

III. INTERDEPENDENT NETWORKS WITH FEEDBACK-DEPENDENCY- 
LINKS 

As mentioned earlier, the above formalism assumes "no-feedback condition", i.e., if a 
node At in network A depends on a node Bi in network B, node Bi may depend only on 
node Ai and not on any other node in network A. In this section we analyse the case 
of interdependent networks without such a constraint. The dependency links between the 
networks are chosen randomly and thus, dependency chains exist, A4 depends on Bi and Bi 
depends on Aj and so on. . 
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In this case Eq. ([TJ becomes [8] 



x -p(l - q) = pqPoo(x). (15) 

where Poo{x) is the relative size of the giant component compared to the initial size of the 
network after random removal of a fraction 1 — x of the nodes. Equation ([2]), the condition 
for first-order transition becomes 

\=pqP' 00 [x). (16) 

Figure compares the case of interdependent lattices with feedback-dependency-links to 
the case of no-feedback condition, which is discussed above in section [H It is seen that the 
feedback case is more vulnerable compared to the no feedback case. The critical threshold 
and the size of the giant component at criticality are both larger compared to the no- 
feedback rule. This can be understood as follows. For small values of q, the probability to 
find a dependency chain of length I exponentially decreases, p(l) ~ q l and the dependency 
chains have a characteristic length, thus, the results are similar to the case with no-feedback 
condition. The length of the dependency chains diverges for q — > 1, where dependency 
chains of the size of the system appear and the system become extremely unstable, p% — 1. 




FIG. 5: The effect of feedback dependency links, a, The critical threshold pc for the 
case of interdependent lattices with feedback-dependency-links (squares) and without feedback- 
dependency-links (circles). The symbols represent results of simulations and the solid lines are 
the solutions of Eqs. ([!]) and ([2]) for the case of no-feedback condition, and of Eqs. (I15p and 
(|16|) for the case of feedback condition, b, The size of the collapse in interdependent lattices with 
feedback-dependency-links (squares) compared to the case of no-feedback condition (circles). 
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Similarly, the solution of a starlike network-of-networks, given before by Eqs. ([6]) and 
((7|), becomes for the case of feedback dependency 



xi = plqPoc^-q + lT- 1 (17) 
x 2 = pqPoo{xi) +p(l - q) (18) 

where the giant component of the root lattice is Poo(^i) and the sizes of the giant components 
of the n — 1 peripheral lattices is -Poo^)- 
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